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Abstract — This work investigates the maximum broadcast 
throughput and its achievability in multi-hop wireless networks 
with half-duplex node constraint. We allow the use of physical- 
layer network coding (PNC). Although the use of PNC for unicast 
has been extensively studied, to our best knowledge, there has 
been little, if any, prior work on PNC for broadcast. Our specific 
results are as follows: 1) For single-source broadcast, the theo- 
retical throughput upper bound isn/(n + l), where n is the 
min-cut size of the network. 2) In general, the throughput upper 
bound is not always achievable. 3) For grid networks, by adopting 
(n+l)-color partitioning and using PNC, the throughput upper 
bound n/(n + 1) is achievable. Our work can be considered as 
an attempt to understand the relationship between max-flow and 
min-cut in half-duplex broadcast networks with cycles (there has 
been prior work on networks with cycles, but not half-duplex 
broadcast networks). 

Index Terms — Broadcast, wireless networks, throughput, 
physical-layer network coding. 

I. Introduction 

This work investigates the maximum broadcast throughput 
and its achievability in multi-hop wireless networks with half- 
duplex node constraint. It is known that in a single- source 
unicast network, the max-flow from the source to the sink 
equals to the min-cut between them. In a single source multi- 
cast network, the max-flow is also equal to the min-cut with 
the adoption of network coding. However, the two results both 
assume full-duplex links that operate independently without 
mutual interference. Our work can be considered as an attempt 
to understand the relationship between max-flow and min-cut 
in networks with half-duplex links that may interfere with 
each other, such as those in wireless networks in which the 
transmission on one link may interfere with transmissions on 
nearby links, and a node does not transmit and receive at the 
same time. 

We allow the use of physical-layer network coding (PNC) 
[1]. PNC is a technique that makes possible the utilization 
of interfering signals. Interference of signals happens when 
there is more than one stream of electromagnetic wave reach 
a receiving node at a given time; most of today's wireless 
networks are designed to avoid interference and interfered 
signals are usually discarded. Rather than shunning interfer- 
ence, embracing interference and transforming it to network 
coding operations at the physical layer is the main thrust of 
PNC. Specifically, when multiple transmitters transmit simul- 
taneously, what is received at the receiver is a superposition of 



the signals. This superposition can be transformed to network- 
coded signals for use in a system that exploits network coding. 

Although the use of PNC for unicast has been extensively 
studied, to our best knowledge, there has been little, if any, 
prior work on PNC for broadcast. Our specific results are as 
follows: 

1) For single-source broadcast, the theoretical throughput 
upper bound is n/(n + 1), where n is the min-cut size 
of the network. 

2) In general, the throughput upper bound is not always 
achievable. 

3) For grid networks, by adopting (n + l)-color partitioning 
and using PNC, the throughput upper bound n/(n + 1) 
is achievable. 

The remainder of this paper begins with the introduction of 
some related works in Section II. In Section III, we analyze 
the broadcast throughput upper bound in a half-duplex wireless 
network. Section IV proposes a scheme to achieve this upper 
bound for grid networks. Finally, section V concludes the 
results in this work. 

II. Related Works 

In graph theory [2], the max-flow min-cut theorem speci- 
fies that the maximum throughput in a single- source unicast 
network is equal to the min-cut. Network coding [3] provides 
a solution to achieve the upperbound min-cut throughput in 
a single-source multicast network. Ref. [4] proposed a model 
for studying the throughput per node of wireless networks. 
Ref. [5] analyzed the throughput of IEEE 802.11 based ad 
hoc networks. PNC, first proposed in [1], incorporates signal 
processing techniques to realize network coding operations at 
the physical layer when overlapped signals are simultaneously 
received from multiple transmitters. It is a foundation of our 
investigation here. Most existing works on PNC focus on the 
unicast scenario. For example, [1] studied the unicast in a two- 
way -relay channel, line networks and 2D grid networks; [6] 
and [7] study the unicast in general networks by designing 
distributed MAC protocols. As far as we know, there has 
been little, if any, prior work on broadcast with physical-layer 
network coding. 

III. PNC Broadcast Throughput Analysis 

We consider the one-source broadcast scenario in which the 
packets from one source X need to reach all nodes in a packet- 



based wireless network. Information of these packets needs to 
be relayed to nodes that are not within the transmission range 
of X by other nodes. 

We represent a packet-based wireless network by an undi- 
rected loopless graph Q = (V, £), where V is the set of nodes 
and £ is the set of links. There is a link {k^k^} G £ between 
two nodes /c, /c' G V if and only if nodes k and k^ are within 
the direct transmission range of each other, i.e., they are within 
each other's one-hop neighborhood. We assume that each link 
has equal capacity. 

Packets generated and transmitted from X are refered 
to as "native packets". We assume equal-sized packets and 
synchronized time-slotted transmissions, in which all nodes 
are scheduled to transmit at the beginning of a time slot. A 
time slot is the duration of one packet. A packet x generated 
and transmitted from X is an element of GF(2*), where s 
is the number of bits in the packet. In other words, x is a 
length- s vector of bits. 

Let JV{k) be the set of adjacent nodes of node k e V. 
Specifically, k' G JV{k) if and only if there is a link {k^k'} G 
£. Each node in our network is half-duplex, i.e., it can either 
send or receive at a given time, but not both. Whenever a node 
k transmits, the same information reaches all its adjacent nodes 
JV{k). We assume that there is no interference from nodes that 
are two or more hops away and there is no transmission loss 
or error. ^ 

We allow the use of PNC reception as defined below. 

Definition 1 (PNC Reception). At the physical layer, when a 
node k receives the superposition of multiple signals contain- 
ing packets yk' G GF(2^) transmitted by several neighbors 
k' G N{k), the node is able to map the composite signal to 
^ = ®k'eN{k) Vk'- We call z a PNC packet. 

Readers who are interested in how PNC reception can be 
realized (with and without channel coding) are referred to [1], 
[8], [9] for details. We remark that yk' transmitted by each 
neighbor k' G M{k) can be either a native packet or a network- 
coded packet. The definition of equation/packet is as follows. 

Definition 2 (Equation/Packet). An equation or a packet is 
expressed as 

z ^ ajXj aixi a2X2 (1) 

3 

where aj.,Xj G GF{2^). It is a linear combination of one or 
more native packets, where aj are the coefficients and Xj are 
the native packets from X. Each of ajXj in (1) is an 5-bit 
vector. Similarly, the packet z in (1) is also a 5-bit vector. If 
there are more than one native packet combined in z, then we 
call it a network-coded packet. 

^This assumption is made to simplify the analysis. In practice, we could 
either employ a forward error control (FEC) scheme or an automatic-repeat- 
request (ARQ) scheme to ensure reliable communication. Both incur some 
overhead in the amount of data to be transmitted. For PNC broadcast, FEC 
is perhaps simpler in that acknowledgements from multiple receivers may 
complicate the ARQ design. 



In accordance with the addition and multiplication opera- 
tions in GF{2^), the addition in (1) is the bit-wise XOR 
over ajXj for different j, and the multiplication of aj and 
Xj for each j in (1) is the multiplication of their polynomial 
representations modulo an irreducible reducing polynomial. 

The terms "packet" and "equation" will be used inter- 
changeably in this paper. A node that has received a packet 
(i.e., z) also has acquired the associated equation, assuming 
the coefficients aj and the identities of the native packets Xj 
(i.e., the indexes of the native packets) are known. Such index 
information can be encoded into packet headers. Conceptually, 
a node will be able to decode D native packets if it has 
D linearly independent equations (native or network-coded 
packets) which contains only the D native packets in their 
summands. 

With reference to Definition 1 and Definition 2, for a 
node k, suppose that a subset of neighbors J\f'{k) C J\f{k) 
transmit simultaneously, and neighbor k' G M'{k) transmits 
Vk' = (^iuf^iuf Then, PNC reception allows node k to 
obtain the following equation. 

^= 2/fe'= l^a,^,x,J (2) 

k'eN'{k) k'eM'{k) \ jy j 

Note that node k' that transmits y^i = 0^^, Xj^, 
may perform upper-layer network coding to obtain y^' = 

©• aj^.Xj,, from the data that it has received so far. That is, 
J/g/ Jk' J k' 

in (2) above, yk' is a packet or network-coded packet generated 
by upper-layer network coding at a transmitting node; whereas 
z is a physical-layer network-coded packet generated at the 
receiving node based on the simultaneous signals received 
from multiple transmitting source. In this paper, we will be 
using PNC as well as upper-layer network coding to enable 
efficient broadcasting. 

Definition 3 (Trivial Network). Consider a network Q = 
(VjS). If all nodes in the network are neighbors of X (i.e., 
A/'(X) = V\{X}), then the network is called a trivial network. 

In a trivial network, all non- source nodes can receive 
directly from X. Relaying information from a node to the 
other is thus unnecessary. There is little value to consider 
trivial networks in our analysis, as the optimal strategy of 
transmission/reception is trivial, i.e., X transmits all the time, 
and the normalized broadcast throughput is 1. In this paper, 
we are only interested in non-trivial networks. 

Consider a cut C = (Vi, V2) in a network Q = (V, S) that 
partitions the nodes V in into two subsets Viand V2 . Let Vi 
be the subset that contains X. Let V[ C Vi be nodes in Vi that 
has neighbors in V2, and V2 ^ V2 to be nodes in V2 that has 
neighbors in Vi. That is, the nodes in Vi\Vi and the nodes in 
V2\V2 are not connected. The information from X has to go 
through some nodes in V[ in order to reach nodes in V2. . 

Definition 4 (Vertex-Cut Size). Consider a cut C. Let 
Vi,V2,V{ and V2 be define as above. The vertex-cut size of 
Cis 




Fig. 1: Partitioning with a qualified cut 



Note the difference between the definition of the traditional 
cut size and the above vertex-cut size. The traditional cut size 
is defined to be the number of edges from nodes in V[ to nodes 
in V2 . The motivation for the above vertex-cut size is due to the 
wireless node constraint we assume: specifically, a node cannot 
transmit different information on different links incident to it; 
when it transmits, it broadcasts the same information on all 
these links. Thus, the vertex-cut size better characterizes the 
maximum flow that can go from Vi to V2. 

Definition 5 (Quafified Cut). Consider a cut C = ( Vi , V2 ) . 

Let Vi , V2 , V( and V2 be defined as above. C is said to be a 
qualified cut if and only if X ^ V{ . 

Fig. 1 shows the partitioning with a qualified cut. A qualified 
cut ensures the adjacent nodes of X are in the same sub- 
network, Vi, as X. A qualified cut does not exist in a trivial 
network, and can always be found in a non-trivial networks. 
Intuitively, the broadcast throughput from X to all other nodes 
is limited by the need to relay information to nodes that 
are not direct neighbors of source X. Thus, the throughput 
limit should be characterized by the qualified cut in that it 
characterizes the "relay capacity" to nodes that are two or 
more hops away from X. 

In this paper, we consider the problem of source X broad- 
casting D native packets to all other nodes in the network 
for large D. Let Wd be the number of time slots (or time 
window) before all nodes acquire all the D native packets. The 
D native packets can be obtained if a node has received D 
linearly independent equations relating the D native packets. 

Definition 6 (Throughput). The broadcast throughput p of 
source X is limi:)^oo D/Wd- 

Each qualified cut has an associated vertex-cut size. The 
qualified cuts with the minimum vertex-cut size in the network 
are called the minimum qualified cuts. The following theorem 
gives an upper bound on the achievable broadcast throughput: 

Theorem 1. Consider a non-trivial network Q whose minimum 
qualified cuts have vertex-cut size n. Then p < n/{n -\- 1). 

Proof of Theorem 1: With respect to Definition 5, let C 
be a minimum qualified cut with vertex-cut size n. Suppose 
the D packets to be broadcast by X are {x{0), ...,x{D — 1)}. 
There are two ways to deliver information related to these D 



native packets from nodes in V[ to nodes in V2: 

1) deliver the information directly in the original form of the 
D native packets; or 

2) deliver at least D linearly independent equations (native or 
network-coded packets) with x(0), x{D — 1) as unknowns. 
Either way, there must be at least D transmissions (taking up 
D time slots) from nodes in V( to nodes in V2. 

Recall that Wd is the total number of time slots needed for 
all nodes to receive D linearly independent equations relating 
the D native packets. We label the nodes in V( by 1, 2, n 
(see Fig. 1). Let Ti,T2,...,Tn be the numbers of packets 
transmitted by nodes 1,2, ...,n, respectively, by the end of 
Wd time slots. Let i?2, be the numbers of packets 

received by nodes 1, 2, n, respectively, by the end of Wd 
time slots. Then with the half-duplex constraint (nodes cannot 
transmit and receive at the same time). 



Wd =Ti-^ Ri for all nodes i,l < i < n. 



(3) 



Furthermore, since each node i^l < i < n, must be able to 
decode the D packets from the receptions (i.e., it must have 
received at least D linearly independent equations), we have 



> D 



(4) 



for all 1 < i < n. In addition, at least D linearly independent 
equations must be delivered to nodes in V2 from nodes in V[ , 
meaning 

Ti+T2 + ...+T^>I). (5) 



Thus, the network throughput is 
D , D 



lim 



lim 



D 



..... ..... < lim 

D^oo Wd d^oo Ti + Ri D^oo Ti^ D 



(6) 



Note that Ti,T2,...,T^ and i?2, (i-^-. how many 

times each node transmits and how many times each node 
receives) depend on the detailed scheduling and relaying 
scheme. Here, we are interested in an upper bound that is valid 
for all schemes, including the optimal scheduling scheme. To 
do so, we can solve the following optimization problem: 

max min lim ^ ^ _ (7) 



l<i<n D^oo Ti -\- D 



such that 



T1+T2 



> D. 



As ^ 00, the solution to the above problem is 

Ti = — ,V1 < z < n. 
n 

The upper bound of p is therefore 

p < lim = 



(8) 



(9) 



The upper bound of broadcast throughput is not always 
achievable. For example, the minimum qualified cut in the 
network shown in Fig. 2 has vertex-cut size 2. However, a 
throughput of 2/3 cannot be achieved. The reader is referred 
to the Appendix for details. 




Fig. 2: Example of a network whose broadcast throughput 
cannot reach the upper bound n/(n + 1) 



Although not always achievable in general, the throughput 
upper bound n/(n + 1) is achievable in grid networks. The 
next section shows how it can be achieved. 

IV. Broadcast Scheme for Grid Networks 

We now argue that throughput of n/(n + 1), where n is 
the min qualified cut size, in grid networks. In grid networks, 
n = 2. To demonstrate its achievability, we have to come 
up with a scheduling scheme specifying for each and every 
node, (i) the time slots during which it transmits; and (ii) the 
information it will transmit during its transmission time slots. 

For the scheduling scheme we will construct, the time slots 
are divided into three disjoint groups (colors) of equal size. 
Each node is assigned one of the three colors. For that the 
source node X, each non-source node can only transmit during 
the time slots of its color. 

We first review the concept of Hamiltonian cycle and present 
a scheme to color the nodes in a grid network based on an 
embedded Hamiltonian cycle. The coloring in turn helps us 
construct the transmission schedules of the nodes. 

A. Coloring of networks by constructing Hamiltonian cycles 

A Hamiltonian cycle is a path that visits each node in a 
graph exactly once and ends at its starting point. First, for any 
M X N grid graph with at least one of M^N being even, 
there is a Hamiltonian cycle in the graph [2]. For example, 
we can construct a comb- shaped Hamiltonian cycle as shown 
in Fig. 3a. Starting from a neighbor of X on the cycle, we 
number along the path by l-2-3-...-(MiV-l). Depicted in Fig. 
3b is an example of this numbering scheme in a 6 x 5 grid 
with X = (1,1). 

Next, if M and N are both odd, then it is not possible 
to construct a Hamiltonian cycle in the grid [2]. Instead, we 
construct a pseudo Hamiltonian cycle called the "Split-Merge 
Hamiltonian cycle", as illustrated in Fig. 4a and 4b. As shown 
in Fig. 4b, after visiting node 5, the visits split into two paths. 
Nodes 6 and 6* are visited in parallel next; node 7 and 7* after 
that; an then node 8 and 8*; finally the parallel visits merge 
back to node 9. By splitting and merging as such, a result is 
the insertion of a new row (as indicated red in Fig. 4b) into 
an (M — 1) X A/" network which already has a Hamiltonian 
cycle constructed (because M — 1 is even, we can construct 
such a cycle). 
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(a) A Hamiltonian cycle in an M x 
network (M is even). 

Fig. 3: M xN networks with at least one of M, being even 
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(b) Numbering of nodes in a 6 x 
5 network after renumbering. 
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(a) A Hamiltonian cycle, part of 
whom is a parallel path, in an 
M X N network (M and are 
both odd). 



(b) Numbering of nodes in a 7 x 
5 network. 



Fig. 4: M X networks with both M, N being odd 



Hamiltonian Node Coloring: Now that we have a number- 
ing scheme for all grid networks, we partition the non-source 
nodes into nodes of three distinct colors c G {0,1,2}. We 
assign a node with number k the color c = k mod 3. We also 
let k* mod 3 := k mod 3, i.e., let the node k* in the split 
paths have the same color as node k. 

The motivation of the coloring scheme based on the Hamil- 
tonian cycle is as follows. A node k that is two or more hops 
away from X can only receive the broadcast information from 
X through its neighbors JV{k). In our time-slotted scheme 
to be described, the nodes with the same color transmit at 
the same time. The Hamiltonian coloring scheme ensures that 
each node k has at least two neighbors assigned with the two 
colors different from the color of node k. These two other 
colors correspond to the time slots in which node k receives. 
Thus, each node k receives in at least two time slots out of 
every three time slots. Section IV-B specifies this transmission 
scheme more exactly. 



B. Ternary transmission schedule 

Definition 7 (Round). A round t^t 
three time slots {3t, 3t + 1, 3t + 2}. 



0, 1, 2, ... is a set of 



Our basic idea is to let nodes with color c transmit in 
time slots 3t + c, t = 0, 1, 2... As a consequence, every node 
transmits once and receives twice in every round. 

If we could ensure that in each round, every node receives 
two packets that contain new information, then the broadcast 
throughput would be 2/3, which is the upper bound given by 
Theorem 1 . Toward that end, we propose a schedule with the 
following rules (also summarized in Table I): 

1) Let the sequence of native packets to be broadcast by 
X be {xo(t)}^=o,i,2,.. and {xi(t)}^=o,i,2,... 

2) In time slot 3t, X transmits xo(t); in time slot 3t + 1, X 
transmits xi{t)\ in time slot 3t + 2, X transmits x^it) 

3) This rule describes the behavior of a node that is not 
adjacent to X\ 

In time slot 3t + c, c G {0, 1,2}, node k ^ M{X), 
with color c transmits 

y\t)=au{z^^{t)®z\{t)), 
where e GF(2**)\{0} is a coefficient, and 
zUt)=rUt-l)- aiy\t-2) 

l:leJ\f(k),l mod 3=(c-l) mod 3, 
1^1 or MN-1 



z1{t)=r1{t-l)- 



aiy''{t-2). 



l:leAf{k),l mod 3=(c+l) mod 3, 
1^1 or MN-1 



In the above, Tq (t — 1), (t — 1) are the two packets 
node k receives from its neighbors in the previous round 
in time slots 3t-3 + ((c- 1) mod 3) and 3t-3 + ((c + 
1) mod 3). (note that these two time slots are the two 
time slots with different colors than c). In the above 
equation, by subtracting aiy^{t — 2), we subtract out 
the information node k sent to an adjacent node / two 
rounds earlier, who now sends back y^{t — 2) multiplied 
by its coefficient ai . This prevents the same information 
from looping back and forth between the node and its 
adjacent nodes. 
4) This rule describes the behavior of node 1 and node 
(MN-l): 

Node 1 and node (MN-l) only transmit native pack- 
ets times a coefficient. Specifically, node 1 transmits 
xo{t - 1) in time slot 3t + 1. Node (MN-1) transmits 
xi (t - 1) in time slot 3t + c, where c = (MN - 1) mod 3 
is its color. We refer to these two nodes as ''virtual 
sources''. Node 1 and node (MN-l) are responsible 
for forwarding {xo(t)}t^o,i,2,.. and {xi(t)}t^o,i,2,.., 
respectively. Henceforth, let node 1 be denoted by Xq 
and node (MN-1) be denoted by Xi. We show below 
that being adjacent to X, they can both derive xo{t — 1) 
and xi{t — 1) by the end of round t — 1. 
In the three time slots in round t — 1, source X transmits 
xo(t — 1), xi{t—l), and xo(t — l)0xi(t — 1) respectively. 
Since each adjacent node of X is colored with one color 
only, it is in receiving mode in two of the three time 
slots. Both xo{t — 1) and xi (t — 1) can be derived based 
on the receptions in these two time slots. 



22 


J.3 


^24 


q25 


26 






q29_ 


^28 


J27 


I20 


J 


^ 


^ 


4 




8 


^ 


^^6 


J5 


I18 


I9 




11 


12 


U7 


3I6 


15 


14 


Jl3 



0/1 






07 

^Z / 


08 
Zo 


I23 




^^1 , 


^30 , 


^29 


I22 ^ 


^ , 


^ , 


-^3 


4 


I2I 


.8 


7 


6 , 


,5 


20 , 


9 


8* 




6* 


19 ^ 


,10 . 


, 


J2 


13 





07 


J6 


J5 


44 



(a) Equivalent of a 6 x 5 network (b) Equivalent of a 7 x 5 network 

Fig. 5: Equivalent networks after removing X 
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Fig. 6: A 2 X 3 network example 



5) This rule describes the behavior of adjacent nodes of X 
who are not the the virtual sources, Xq and Xi. By this 
rule, we ensure that only the two virtual sources can 
send out the newest native packets of X: 
An adjacent node k of X, who is not Xq or Xi, can 
also derive xo{t — l) and xi(t — 1) by the end of round 
t — 1. In round t, node k transmits 

where a/^ G GF(2*)\{0} is a coefficient, and 

^oW = ^oW-4(i-l) 
znt)=z'[{t)-x[{t-l). 

In the above, ZQ(t),zf(t) are defined as in rule 3, and 
Xo(t — 1) and x[(t — 1) are the two packets sent by X 
and received by this node in round t — 1 . 
For example, if anode k has color 0, then XQ{t — 1) 
and x[{t — 1) are the two packets sent by X in time 
slots 3t — 1 and 3t — 2, respectively; i.e., XQ(t — 1) = 
xo{t - 1) © xi(t - 1) and x[{t - 1) = xi(t - 1). 
Transformation to Two-Source Broadcast Problem: With 
the above rules, only the virtual sources Xq and Xi can 
be considered as the origins of the newest information. The 
single-source networks in Fig. 3b and Fig. 4b can then be 
transformed to two- source networks in Fig. 5a and Fig. 5b. In 
Fig. 5a, nodes 1 and 29 are Xq and Xi, respectively; in Fig. 
5b, nodes 1 and 31 are Xq and Xi, respectively. 

Example 

We first illustrate what happens when applying this schedule 
to a simple 2x3 grid network. The source is located at 
(0,0). Fig. 6a shows the numbering of this network and Fig. 
6b shows the corresponding coloring. In this example, we 
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TABLE I: Ternary transmission schedule 



set all coefficients to 1. The transmissions of our ternary 
schedule are shown in Table 11. Although node 1 and node 4 
can overhear each other, they are of the same color and thus 
transmit at the same time. Hence they will not interfere each 
other. Note that in time slot 9, node 3 transmits xo(0) 0xi(O). 
Upon receiving xo{0) xi{0), node 2 can decode xi{0) 
because it already has xo(0); node 4 can decode xo{0) because 
it has xi{0). From round 3 onwards, every non- source node 
receives sufficient information for it to derive two new native 
packets in each round. Thus, the throughput is 2/3. 

It can be observed that the two virtual sources node 1 and 
node 5 always have (0) , . . . , (t) and xi (0) , . . . , xi (t) by the 
end of round t. Therefore they can always derive xo{t -\- 1) 
and + 1) by the end of round t -\- 1 even if they receive 
PNC packets during round t + 1, because they know all but 
one of the unknowns (native packets). 

In this example each node can obtain two native packets 
in a round. In a general network where there is interference 
among nodes, obtaining native packets cannot be guaranteed. 
However, as will be shown, we could ensure every non- source 
node still obtains two linearly independent equations in each 
round. 

In a general grid network, depending on its position in 
the grid, a node can have up to four neighbors. With the 
Hamiltonian Node Coloring, it is possible for a node to 
have two or three neighbors of the same color, only one of 
which is the adjacent node on the Hamiltonian cycle (see 
Section IV- A on Hamiltonian Node Coloring). When multiple 
neighbors of the same color transmit simultaneously, the node 
receive a PNC packet, for which the XOR of the simultaneous 
transmissions of the neighbors of the same color is received. 
For example, in Fig. 4b, node 9 (color-0) has four neighbors: 
node 8 (color-2), node 8* (color-2), node 10 (color- 1) and 
node 20 (color-2). As a consequence, in each round node 9 
receives from three nodes simultaneously in the color-2 time 
slot, which yields a PNC packet; it receives from only one 
neighbor in the color- 1 time slot. For both packets, we need 
to make sure: 

1) the packet received is linearly independent with all 
packets previously received; 

2) the packet received is not null. 

C. Random Coefficients 

In the previous simple 2x3 example, the coefficients ak 
for all nodes k, k ^ 1 or MN — 1, were set to 1. In a general 
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TABLE II: Ternary transmission schedule for a 2 x 3 network, 
"s", "r" and "d" indicate "send", "receive" and "derive", 
respectively. 



grid, this scheme may not work. Henceforth, we consider 
a time- varying random coefficient scheme. Specifically, the 
coefficient ak{t) of node k in round t is chosen uniformly 
randomly from the non-zero element of GF{2^), and ak{t) 
for different k and t are i.i.d. The sequences of packets 
xo(t)^^Q^2 ^i(^)t=oi2 transmitted by the virtual 

sources Xq and Xi remain the same, and their coefficients 
can be considered as 1 throughout the process. We show that 



with high probability, things will work out so that the grid 
network will have throughput of 2/3. 

In grid networks, the newest information (i.e., signals em- 
bedded with the latest native packets) could come through a 
shortest path from the the virtual sources Xq and Xi rather 
than through the sequence of nodes along the Hamiltonian 
cycle. 

Definition 8 (Siiortest Patii). A shortest path from a virtual 
source to a node is a shortest sequence of adjacent nodes 
leading from the virtual source to the node, in which no two 
adjacent nodes in the sequence are of the same color. 

Note that in general there could be multiple shortest paths 
of the same length leading from a virtual source to a node, and 
some of them may share some common intermediate nodes. 

Definition 9 (Hamiltonian Sub-paths). The paths from the 
virtual sources Xq (numbered 1) and Xi (numbered MN — 1) 
to a node k along the Hamiltonian cycle, i.e., po = 1 — 2 — 
... - k and pi = {MN - 1) - {MN - 2) - ... - k, are called 
Hamiltonian sub-paths of node k. 

Definition 10 (Coefficient Product of a Path). The coefficient 
product of p = ki — k2 — ... — ki — k, 3. path from node ki to 
node k, in round t is 

gp{t) =akAt-{i-i))akAt-{i-^)) 

...ak,_At-l)akM (10) 
I 

= l[(^kAt-{l-r)). (11) 

r=l 

where {t-{l- 1)), {t-{l-2), aki_^{t - 1), Offc^ {t) 
are the coefficients of nodes /ci, /c2, /c^-i, /c^ in rounds t — 
(/ — 1), t — (/ — 2), t — 1, respectively. 

If a native packet x goes through a path p, then its coefficient 
when it arrives at the last node will be gp{t), the coefficient 
product of this path. If a native packet begins its journey from 
a node j with an initial coefficient aj and arrives at a node k 
in the same time slot in round t via multiple paths Pi,P2, 
then the coefficient of the native packet x when it is received at 
node k is the sum of coefficient products otjQp^ {t)^ajgp^ (t)0 
. . . Note that these multiple paths must have the same number 
of hops since different copies of packet x are supposed to 
arrive at node k together in the same time slot. 

Lemma 1. The coefficient product gp{t) of a path p = ki — 
k2 — ... — ki — k is non-zero for t > I — 1, and is uniformly 
randomly distributed over in GF{2^)\{0}. 

Proof of Lemma 1: Since a/ci(t — {I — l)),a/c2(t — {I — 
2), Q^/e^_i — 1), Qf/c^ (t) are i.i.d. uniformly distributed over 

GF(2^)\{0} when t>l-h gp{t) = Ul=i o^kA^ - - r)) 
is non-zero and is uniformly distributed over GF{2^)\{0}. ■ 

Example of received packets 

For example in Fig. 4b, in time slot 3t, node 10 receives 
from node 9. The shortest path from Xq to node 10 is = 1 — 
8—9—10. However, this is the shortest path for {xo(t)}^^o,i,..., 



because node 1 only transmits {Q^i(t)xo(t)}^=o,i,...- Shortest 
paths from Xi to node 10 for {xi(t)}t=o,i,... are pi = 31 — 
2 - 1 - 8 - 9 - 10, P2 = 31 - 2 - 7 - 8 - 9 - 10 and 
Ps = 31-2-7-8* -9- 10. A native packet in {xo{t)}t=o,i,... 
will be multiplied by the coefficient of the sender when being 
sent. The packet node 10 receives from node 9 in time slot 3t 
is 

^o'W =gpo{t)xo{t - 3) 

e {gp^ {t) © gp, {t) © gp^ {t))xi {t - 5) 

©/(xo(t-4),xo(t-5),...,xo(0)) 
©/(xi(t-6),xi(t-7),...,xi(0)) (12) 
=ai{t - 2)as{t - l)aAt)xAt - 3) 

© (a3i(t - 4)a2(t - 3)ai{t - 2)aAt - l)aAt) 
© a^i{t - 4)a2{t - 3)a7{t - 2)as{t - l)ag{t) 
© asi{t - 4)a2(t - 3)a7{t - 2)as* {t - l)ag{t)) 
xi{t - 5) 

©/(xo(t-4),xo(t-5),...,xo(0)) 
©/(xi(t-6),xi(t-7),...,xi(0)) (13) 

where /(•) is a linear combination of the arguments. 

The time indexes of the newest native packets in Eqn. 
(13) escalate over time. Therefore each new rQ^{t) is linearly 
independent with all rl^{t'),t' < t. 

rl^{t), the packet node 10 receives in time slot 3t-\-2 will 
be a PNC packet from nodes 11 and 17. A shortest path in this 
time slot for {xo(t)}t=o,i,... is 1-8-21-20-19-18-17-10; 
and a shortest path for {xi(t)}^=o,i,2,... is 31-2-7-8*-7*- 
12 — 11 — 10. Therefore the newest native packets in it would 
be xo{t — 7) and xi{t — 7), multiplied by some ®p^s 9pi^)' 
where S is the set of all shortest paths and p is a shortest path 
for the native packet. 

Packets received by a general node 

Consider a general node k^k ^ 1 or MN — 1, in the net- 
work. Focus on one of its two receiving time slots in round t. 
Let 5*0 , Si be the sets of shortest paths for {xo(t)}t^o,i,... and 
{^i(^)}^=o,i,.... respectively, in this time slot; and Sq^^S^^ 
be the sets of paths that are q hops longer than and Si for 
{^o(^)}t=o,i,... and {xi(t)}t=o,i,..., respectively. Depending 
on q, each of 5*^^ and S^^ may or may not be empty. For 
example when q < 0, both Sq^ and 5*^^ will be empty and 
thus a^+q {t) = and a^+q {t) = 0. 

In round t, node k receives 

^o(^) =aso{t)xo{t - io) © ag+i{t)xo{t - io - I) 
© ... ©a^+t-io {t)xo{0) 
© asAt)xi{t - ii) © ag+i{t)xi{t - n - 1) 
© ... ©a^+t-.i(t)xi(0), (14) 

where io and ii are the lengths of paths in 5*0 and Si, 
respectively; and 

^So{t) = gp{t). as,{t) = gp{t), (15) 

P^Sq p^Si 



(16) where 



Grouping the packets received in this colored slot in all rounds, 
we have a linear equation system {^o(^)}^>min(io,n)- 

Now consider the other receiving time slot of node k A{D — 1) 
in round t. Let To,Ti be the sets of shortest paths for 
{xo(t)}t=o,i,... and {xi{t)}t=o^i^..., respectively, in this time 
slot; and Tq^^T^^ be the sets of paths that are q hops 
longer than Tq and Ti for {xo(t)}^^o,i,... and {xi(t)}^^o,i,..., 
respectively. Node k receives 



zi{t) =aTo{t)xo{t - jo) © %+i {t)xo{t - jo - 1) 
© ... ©a^+t-io {t)xo{0) 

^ 

© ttTi {t)xi{t - ji) © aj.+i{t)xi{t — ji — 1) 

© ... ©a^+t-n(t)xi(0), (17) 
^ 1 

where jo and ji are the lengths of paths in Tq and Ti, 
respectively; and 



uo(A:-l) ui(A:-l) 

uo{D^k-2) ui{D^k-2) 

vo(fe'-l) vi(fc'-l) 

_ vo(I^ -^k' -2) vi{D -^k' -2) 



(23) 



- xo(0) - 










xo(i:>-i) 

XI (0) 


z(i^-l) = 


zo{D + k- 
zi(k' -1^ 


2) 

) 


• (24) 


_ xi(D-l) . 




. zi{D + k' - 


2) - 





In ^4(1^ — 1), uo(t), ui(t), vo(t), vi(t) are defined as follows: 



pen 



gp{t). 

peTi 



(18) 



a^+q (t) 



gp{t), a^^,(t)= ^^(t). (19) 



uo(t) - 



Grouping the packets received in this colored slot in all rounds, 
we have another linear equation system{zi(t)}^>min(jo,ji)- 

For simplicity, we assume the number of native pack- 
ets from the source X is an even number 2D, where 
D is an integer. Hence there are D native packets from 
each of the virtual sources Xo and Xi. As a consequence, 
{^o(^)}^=o,i,... = {^o(^)}^=o,i,...,D-i and {xi(t)}t=o,i,... = 
{^i(^)}^=o,i,...,D-i- After all native packets have been sent by 
the source, we allow MN more time slots for them to circulate 
in the network. During these MN time slots, the source and 
the virtual sources can be considered as transmitting null 
packets, which do not increase the number of unknowns in 
the network. For node k,l < k < MN — 1, let 



CLa + t-iQ (t) 

a +t-i-iQ (t) 

a^+o(t) 




ui(t) = 



a^+t-ii (t) 
a^+t-i-ii (t) 

a^+o(t) 



(25) 



vo(t) 



a +t-jQ (t) 
a +t-i-jo (t) 



ttrp + O (t) 

°0 







VI (t) 



a +t-3i (t) 
^ 1 



ttrp + O (t) 







(26) 



k' = MN - k. 



(20) 



We select {zo{to)}k-i<to<D+k-2, a subset of 
{^o(^)}t>min(io,ii)' and {zi{ti)}k'-i<ti<D+k'-2, a subset 
of {^i(t)}t>min(jo,ji)' group them together as a single 
linear equation system 



{^o(^o), Zi{ti)}k-l<to 



<D+k-2,k' -l<ti<D+k' -2- 



(21) 



If all equations in (21) are linearly independent, then we 
can derive all 2D native packets from (21). For this purpose, 
we represent (21) in matrix form as follows: 



A(D- l)x = z(D- 1), 



(22) 



Note that by definition 5'+^ 51+^ To+^ are the shortest 
paths 5'o, To, Ti, and ^o,ii, jo, ji are the lengths of paths 
in 6*0, To, Ti, respectively (defined in (14) and (17)). There 
are max(0, T) — t + io — l),max(0, T) — t + ii — 1), max(0, T) — 
^ + jo — 1) and max(0, D — t ^ ji — 1) trailing entries in 
uo {t) , ui {t) , vo (t) , vi {t) , respectively. 

Before we discuss the relationship between the coefficients, 
we have some definitions. In the following, we will use \p\ to 
represent the length of a path i.e., the number of hops in p. 

Definition 11 (Equal-hop Path Set). An equal-hop path set 
is a set of paths of the same length. For example S = 
{pi,P2, •••}, where \pi\ = \p2\ = is an equal-hop path set. 

Definition 12 (/i-hop Path Set). An equal-hop path set S is 
called an h-hop path set if all paths p e S have the same 
length \p\ = h. 



Lemma 2. Consider two equal-hop path sets S = 
{PiiP2-, '"-tPq] cind S' = {Pi,P2 5 •••5-Pr}- Suppose that there 
is at least one path in S that is not in S' or there is at 
least one path in S' that is not in S. Let d G GF(2^) be 
any given fixed value. Let as(t)^t > maxi<^<g(|p^|), and 
as>{t')-,t' > maxi<^<^(|pj|), be the coefficients of a non-zero 
packet after going through S and S' , respectively. Then asif) 
and as'{t') are independently distributed and the probability 
that as{t) = das'(t') is no more than 1/(2* — 1). 

Proof of Lemma 2: Consider a native packet x with no 
coefficient at the beginning. After going through S and 5*^ the 
coefficients of the packet will be 

pes pes' 

t > maxi<^<q(|pi|) and t' > maxi<i<^(|p-|) ensure that 
gp{t) y^O.peS and gp{t') ^ 0,p G S\ 

Without loss of generality, suppose there is one path p in 
5* that is not in S' . The coefficient product of p, gp{t), is 
independent of any other coefficient product in as{t) and 
as'{t'). Thus, as{t) and as'{t') are independently distributed. 
By Lemma 1, gp{t) is non-zero and is uniformly distributed 
over GF(2*)\{0}. The probability for 



i.e., for 



as{t) = das'{t'), 

gp{t) = das'{t') - as\{p}{t) 

= d^gp(t')® gp{t) 



(27) 



pes' 



pes.p^p 



is 1/(2^-1) when ^0^^^, ©©pG5,p/p W ^ 0' ^nd 
when d@^^g, gp{t') © ®p^s,p^p9pif) = 0- Therefore the 
probability that (27) holds is no more than 1/(2^ — 1). ■ 



Lemma 3. Consider the equal-hop path sets Sq 



+/c — 1 — io 



+k'-l-j 



and 

^ .The probability that a^+k-i-i^ (t) ^ Ofor any t > 
k — 1 is no less than 1 — 1/(2* — 1); the probability that 
a^+k'-i-ji{t') ^ for any t' > k' — 1 is also no less than 
1 - 1/(2* - 1). 

Proof of Lemma 3: Recall that is the length of paths 
in S'o and ji is the length of paths in Ti . The length of paths 
in is — 1, which is the length of the Hamiltonian 

sub-path po = 1 — 2 — ... — /c. 

Case L If 5'+^"^"^° has only one path pq, then by Lemma 1 
a +fc-i-io(t) 7^ with probability 1. 

Case II: If S^"-'-'" has more than one path, we divide it 
into two sets 5' = {po} and 5" = 5'^^"^"'°\{po}. Then 
a +fc-i-io(t) = + a^/. By Lemma 2 the probability that 
o^s = OjS' is no more than 1/(2* — 1) when t > k — 1. Therefore 
a +k-i-io (t) ^ with probability no less than 1 — 1/(2* — 1). 

In conclusion, the probability that a^+k-i-iQ (t) ^ for any 
t > k — lis no less than 1 — 1/(2* — 1) . Similarly we can prove 



the probability that a^+y-i-j-^ (f) ^ for any f > k^ — 1 is 
also no less than 1 - 1/(2* - 1). ■ 



Lemma 4. a^^ (t), a^^ (t) in Eqns. (15)-(16) and 
aTo(^), ciTi {t) in Eqns. (18)-(19) are non-zero with probability 
no less than 1 — 1/(2* — 1). 

Proof of Lemma 4: We induct on the coefficients of 
the newest native packet from {xo(t)}t=o,i,... in the packet 
received during round t in the following discussion. 

When node k has a shortest path of zero hop from the virtual 
source Xq, i.e., node k is the virtual source Xq, the shortest 
path is directly from X to node k. It receives xo{t) within 
the newest packet and the coefficient of xo{t) is 1, which is 
non-zero. 

When node k has a shortest path of one hop from Xq, it 
receives ai{t)xo{t — l) within the newest packet, where ai{t) 
is the coefficient of Xq . Since ai {t) is uniformly distributed 
over GF(2^)\{0}, it is non-zero. Now we assume all shortest 
paths of io hops from Xq have non-zero coefficients for xq (t — 
io)- 

When node k has a shortest path of + 1 hops from X. 
Suppose all its adjacent nodes on the shortest paths from X 
have non-zero coefficients for the newest native packet from 
{xo(t)}t=o,i,.... Consider node an adjacent node of node 
k, that is on some shortest paths from Xq to k. We focus 
on one receiving time slot of node k, and two receiving time 
slots of node /. Let the sets of shortest paths from Xq to / 
for {xo{t)}t=o,i,... in the two receiving time slots of a round 
be 5*0, ^ and ^. Suppose node Vs newest native packets from 
{xo(t)}^=o,i,... within the two packets received are 



ao,i{t)xo{t - io,i) and bo,i{t)xo{t - jo,i), 



where 



^oAt) = gp{t), boAt) = gp{t). 

and io,hjo,i are some non-negative integers with 
min(io,^, jo,0 = ^o- 

The newest native packet from {xo(t)}t^o,i,... within the 
packet sent from node / to node k then will be 



ai{t)[ao^i{t)xo{t - io^i - 1)] if io^i < jo,z 

^i{t)[{ao,i{t) © bo,i{t))xo{t - io^i - 1)] if io,i = jo,i , 
^c^i{t)[bo,i{t)xo{t - jo^i - 1)] if io,i > jo,i 



Because ao,^(t) and bo^i{t) correspond to different sets of 
shortest paths in different time slots of a round, and from 
Lemma 2, ao^i{t) ^ bo^i{t) with probability no less than 1 — 
1/(2* — 1), we conclude ao^i{t) +6o,z(t) / with probability 
no less than 1 — 1/(2* — 1). By assumption, ao,z(t) is non-zero 
if ^0,^ < jo,i, because io,i = min(zo,z, jo,0 = corresponds 
to paths of io hops from Xq; similarly boj{t) is non-zero if 
^0,/ > jo,h CLo,i{t) and boj{t) are both non-zero if io,z = jo,i- 
The probability for the coefficient of xo{t — io-\-l), the newest 



native packet from {xo{t)}t=o,i,... within y^{t), to be zero is 
one of the following: 

1 



FT{ai{t)ao,i{t) = 0) = Pr^W = {ao,i{t))~^) 
FT{ai{t){ao^i{t)ebo^i{t))=0) 
= FT{ai{t) = (ao,K^)© VK^))"') = WT—^ 



2^-1 



FT{ai{t)boAt) = 0) = FT{ai{t) = (boAt))'^) 



1 



These will be the same as the probability that as^it) = 
when there is only one adjacent node that is in the shortest 
paths from X to k and transmitting. 

When there are two adjacent nodes of node k that are in the 
shortest paths from X to k and transmitting together, let these 
nodes be denoted by /q, ^i- Without loss of generality, suppose 
^0,^0 < Jo,Zo and io,^i < jo,h' The received PNC packet will 
be 

z^{t) ={ai^{t)ao,i^{t) © ai^{t)ao,i,{t))x{t - io) 

®f{xi{t-ji),...,xi{0)) 

where io is a non-negative integer. Hence aso{t) = 

C^lo (^)<^0,^o (^) ® (^h {'^)^0M (^)- 

Because azo(^) and ai^{t) are independently uniformly 
distributed over GF(2^)\{0}, the probability for aso{t) = 
is 



Pr(as„(t) = 0) 



1 



When there are three adjacent nodes of node k that are in 
the shortest paths from X to k and transmitting together, let 
these nodes be denoted by lo^h^h- Without loss of generality, 
suppose zo,^o < io,Zo and zo,Zi < jo,h- The received PNC 
packet will be 

© ai^{t)ao,i^{t))x{t - io) 
©/(xo(t-io-l),...,^o(0)) 

e/(^i(^-ji),...,^i(o)) 

where io is a non-negative integer. Hence as'o(^) — 
o^io W«o,^o W ® W<^o,^i it) © Wao,Z2 

Because (^1^(1), ai^{t) and a/2(t) are independently uni- 
formly distributed over GF(2*)\{0}, the probability for 



FT{as,{t) = 0) 



:PrK(t) 



There is no such scenario that there are four adjacent nodes 
transmit together, because there must be one different color 
among them, as^ then will be non-zero with probability no 
less than 1 - 1/(2^ - 1). 



By similar arguments, we can show a^^ , a^o , clti are also 
non-zero with probability no less than 1 — 1/(2^ — 1).. 



Lemma 5. Equations inside (21) are linearly independent with 
probability at least (1 - 2^31) (1 " ^^?r^)^^"^- 

Proof of Lemma 5: Consider the first two equations 
zo{k — 1) and zi{k' — 1) inside (21). The coefficient of xo{0) 
in zo{k — 1) is 



a +k-i-io {k — 1) 



9p{k - 1), 



pes^ 



where 5'^^"^"'° is the {k - l)-hop path set from Xq to node 
k, the paths in which have k — l — io more hops than a shortest 
path, which is io hops in length; and the coefficient of xi{0) 
in zi{k' — 1) is 



ft^+fc'-i-ji {k' — 1) — 



per; 



E 

+k'-i-ji 



9p{k' - 1), 



where 



^ is the {k^ — l)-hop path set from Xi to node 



k, i.e., the paths that have k' — 1 —ji more hops than a shortest 
path, which is ji hops in length. Note that by Lemma 3, with 



probability no less than 1 — 1/(2^ — 1) that 



is 



non-zero; similarly, with probability no less than 1 — 1/(2^ — 1), 

(2 -Lfc/ _i_ j-^ (k^-l) is non-zero. In other words, with probability 

^ 1 

no less than 1 — 1/(2* — 1), xo{0) will appear in zo{k — 1) 
and with probability no less than 1 — 1/(2* — 1), xi{0) will 
appear in zi{k' — 1). 
There are two cases: 

1. xo{0) does not appear in zi{k' — 1) or xi{0) does not 
appear in zo{k' — 1). This means that if xo{0) appears in 
zo{k - 1) and xi(0) appears in zi{k' - 1), then $d G GF{2') 
such that zo{k — 1) = dzi{k' — 1); 

2. both xo(0) and xi(0) appear in zo{k — 1) and zi{k' — 1): 
Note that the Hamiltonian sub-paths (defined in Definition 9) 
are po e 5^^"^"^° and pi e T+^'"^"^'\ However, po ^ 
T(+^"^°,0 <t<k' -h and pi ^ S^*~'\0 < t < k - 1. By 
Lemma 2 the probability that 3d G GF{2^) such that 

a +k-i-io (k — 1) = d - (coefficient of xo(0) in zi(k' — 1)) 

and 

coefficient of xi{0) in zo{k — 1) = d ■ a +^'-1-^1 {k' — 1) 

is no more than 1/(2* — 1). To see the above statement, note 
that we certainly can find a d to satisfy the first equation: 
i.e., d — a^+k-i-io{k — 1) • (coefficient of xo{0) in zi{k^ — 
1))~^. However, given this fixed d, the probability of satisfying 
the second equation is no more than 1/(2* — 1) according to 
Lemma 2. Therefore the probability that zo{k — l) = dzi{k' — 
1), i.e., zo{k — 1) and zi(k' — 1) are linearly dependent, is at 
most 1/(2* — 1). As a consequence, the two equations inside 
each pair of equations in {^0(^0), ^i(^i)}to=/c-i,ti=/c'-i are 
linearly independent with probability at least 1 — 1/(2* — 1). 



We now proceed with an inductive proof on the Hnear 
independence between equations of different pairs. Suppose 
all equations inside 



{zo{to),Zi{ti)}k- 



l<to<k-l-\-A,k' -l<ti<k' -l-\-A^ 



(28) 



for some A, where < A < I) — 1, are linearly independent. 
Similar to the matrix form of (21) in (22), we represent (28) 
in matrix form as follows: 



where 



A(A) 



A(A)x = z(A), 
uo(A;-l) ui(A;-l) 



(29) 



uo(A;-l + A) 
vo(fc'-l) 



ui(fc - 1 + A) 
vi(fc' - 1) 



(30) 



vo(fc'-l + A) vi(A:'-l + A) J 

zo{k-l) 



«o(fc-l + A) 
ziik'-l) 



- xo(0) - 




xo{D-l) 
XI (0) 


z(A) = 


_ xi{D-l) . 





(31) 



in which uo(t), ui(t), vo(t), vi(t) are defined in (25)(26). 
74(A) is a 2(A + 1) x 2D matrix. By assumption all rows 
in A are linearly independent, i.e., A has full rank 2(A + 1). 
The pair of equations just after the set in (28) are 

zo{k + A) =as,{k + A)xo{k + A - zq) 

© a^+i (k + A)xo(A: + A - - 1) 
© ... © a^+fc+A-.o {k + A)xo(O) 
©a5i(A: + A)xi(A: + A-ii) 
©a^+i(A: + A)xi(A: + A-ii -1) 
© ... © a^+fc+A-.i (k + A)xi(O), (32) 

zi{k' + A) =aT^{k' + A)xQ{k' + A - jo) 

© a^+i (k' + A)xo(/c' + A - jo - 1) 
© ... © a^+e+A-io {k' + A)xo(O) 
© ttTi (/c' + A)xi (fe' + A - ji ) 
© a^+i (fe' + A)xi (/c' + A - ji - 1) 

(fc' + A)a;i(0), (33) 

^ 1 

the coefficients of which are represented by vectors 

[ uo(fe + A) ui(fe + A) ] (34) 

and 

[ vo(fe' + A) vi(fe' + A) ], (35) 



respectively. By the same argument 
previous case of A = 0, zo(/c 
zi{k' + A) ( and thus [ uo(A: + A) 
and [ vo(/c' + A) vi(fc' + A) ]) can be 



as 



for the 
+ A) and 
ui(fc + A) ] 
shown to be 



linearly independent with each other. Before deriving the 
probability that equations (28)U{zo(/c + A), zi{k' + A)} are 
linearly independent, we first evaluate the probability that 

3/3 G GF(2^)2(^+i) such that 



/3A(A) = [ uo(fe + A) ui(fe + A) ] . 



(36) 



That is, we consider [ uo(/c + A) ui(/c + A) ] first without 
including [ ^ro{k' + A) vi(/c' + ^) ] • If a vector such as (5 
exists, then [ uo(/c + A) \i\(k + A) ] is not linearly inde- 
pendent with A(A) (i.e., zo(/c + A) is not linearly independent 
with (28)). 

We use A(A)(^^ ,p2,.- ;ci,c2,.- ) denote the submatrix of 
A(A) that consists of the intersection of rows ri,r2,... and 
columns ci,C2,.... Let 

A = A(A)(i^...^2(A+l);l,...,A+l,D+l,...,D+A+l), 

where A is a2(A + l) x2(A + l) square matrix that contains 
the intersection of all rows in A (A) and columns 1 through 
A + 1 and columns D^l through i:> + A + 1 in A(A). 



We use u 



(C1,C2,.--) 



to denote the subvector of u that consists 



of the ci-th, C2-th, ... elements of u. Let 

a = [ uo(/c + A) ui(/c + A) ] . 

(37) 

u is a vector of length 2(A + 1) that contains entries 1 
through A + 1 and entries 1^ + 1 through D + A + 1 in 
[ uo(A: + A) ui(/c + A) ]. If 3/3 G GF(2^)2(^+i) such 
that (36) holds, then 



(3 A 



(38) 



holds, too. Now we evaulate the probability of solving /3 from 
(38). 

The first column of A is 



a +fc-i-io {k — 1) 
a^+k-iQ (k) 



i(l,..,2(A+l);l) 



a^+fc-i+A-io (k 



1 + A) 



^0 



a +y-i+A-jo{k' -1 + A) 

-L n 



Recall that one of the two Hamiltonian sub-paths po is in the 
path set 5^^"^"^° but not in 6'^^"^°,0 < t < A - l,t ^ 
k - 1 or To+^"^°,0 < t < A - 1. By Lemma 2 the first 
entry ft^^+fc— i— iq {k — 1) in 2(A+i);i)^ which is the sole 
coefficient related to po in this column, is independent of any 



other entry in the column. In general, the S-th column of A, 
1< (5 < A + 1 , is 



a ^+k-5-iQ (k — 1) 



(fc-l + A) 

1(1,..,2(A+1);.) = °a ,_,_^^(k'-l) 

^0 

Cirj.+ k' -S+l-jQ ik') 
^0 



a +fc/-5+A-jo (A:' - 1 + A) 

n 



in which by Lemma 2 the (5-th entry a^+fc-i-io(/^ — 2 + (5), 
which is the sole coefficient related to po in this column, 
is independent with any one of the other entries. Thus, we 
conclude that the 5-th entry in column (5, l<5<A + l,of 
A is independent of any other entry in column 5 of A. 
The (A + 2)-th column of A is 

a^+fc-i-ii (k — 1) 
a {k) 



Ail 



(l,..,2(A+l);A+2) 



a„+fc-i+A-ii {k — 
a^+k'-i-i^{k' - 1) 



A) 



n(fc') 



A) 



Recall that Hamiltonian sub-path pi is in the path set 

T+^'"^"^^ but not in 5+'"^^ , < t < A - 1 or T+^"^^ , < 
t < A - l,t 7^ fc' - 1. By Lemma 2 the (A + 2)-th entry 
(/c' — 1) is independent of any other entry. The (5-th 

column of A, A + 2 < 5 < 2(A + 1), is 



^(l,..,2(A+l);5) 



tl^+fc-(5 + A + l-ii — 1) 



Qj j^ + fc-(5 + A + 2- 



at^+fc-(5+2A+i-ii {k — \ -\- A) 

a +fc/-5+A+i-ii (/c' — 1) 
^ 1 

«^ + fc'-<5 + A + 2-ji {k') 



<^^ + fc'-<5 + 2A + l- 



,,(fc'-l + A) 



in which by Lemma 2 the 5-th entry a +fc 

1 



,_,,(A:^ + (5-A- 

3) is independent of any other entry. Thus, we conclude that 
the (5-th entry in column (5, A + 2 < (5 < 2(A + 1), of A 
is independent of any one of the other entries in the same 
column. 



Therefore the diagonal entries of A are independent of 
any other entry in the same column. After triangularization 
by elementary row operations, the resulting diagonal entries 
will be uniformly distributed over GF(2^). If none of these 
diagonal entries are 0, the triangularized matrix of A will have 
full rank. Therefore A has probability (1 - 1/(2^ - 1))(1 - 

1/(2^ - = (1 _ 1/(2^ _ l))2(A+l) + l ^^^^ f^ii 

rank, i.e., to be invertible. 

Since A has probability (1 - 1/(2^ - i))2(a+i)+i 
be invertible, (3 can be solved from (38) by /3 = 
uA-^ with probability (1 - 1/(2^ - i))2(a+i)+i j^^^ 
[ uo(fc + A) ui(/c + A) ] in (36) is fixed to be /3A(A). 



However, uo(/c + A) contains 



,{k + A), the time- 



varying random coefficient of xo(A + 1) in Z{){k-\- /S), which 
is associated with the Hamiltonian sub-path p^, as its (A + 2)- 
th entry. The (A + 2)-th column of A(A) is 



(2(^ + fc-A-2-io ~ 1) 

a^^+fc-A-i-ig (k) 



^(A)(i,..,2(A+l);A+2) = 



a +fc-2-*o (/c - 1 + A) 
ct +y-A-2- jo {k^ — 1) 

^0 

Cirj.+ k' -A-l-jQ {k') 



a +k'-2-jo{k' - 1 + A) 



where every entry is independent of a ^ 



)(A: + A). Hence 



a^+k-i-iQ{k + A) has only 1/(2^ — 1) probability to be the 
value of the (A + 2)-th entry in /3A{A). 

On the other hand, A has probability 1 - (1 - 1/(2^ - 
l))2(A+i)+i j^Q^ invertible. In this case, if u in (38) lies 
in the image of A, then (38) can be satisfied, and there will be 
more than one (3 that can satisfy (38). If u does not lie in the 
image of A, then (38) cannot be satisfied. For an upper bound, 
we can simply assume that whenever A is not invertible, (36) 
can be satisfied. Specifically, overall the probability that 3p G 
GF(2^)2(^+i) such that (36) holds is at most 



1 



2^ - 1 
Since 



(1 



2(A+1) + 1 



+ (1-(1 



2(A+1) + 1 



)• 



1 > (1 _ ^)2(A+iHi > 1 _ 2(A + 1) + 1 

- V 2^-1^ ~ 2^-1 



we have 

( ' 



)(1 



1 



\2(A+1) + 1 



< 



2^-1'' 2^-1' 
1 , (2(A + 1) + 1) 



2(A+1) + 1 



< 



"2^-1 
_2(A + 2) 
'2^-1 
2(1^ + 1) 



1 



Hence 



A' = 



uo(^-l) ui(fe-l) 



uo(fc + A) 
vo(^'-l) 



ui(fc + A) 



_vo(^'-l + A) vi(/c^-l + A) 

has probability at least 1 - 2(Z:> + 1) /(2^ - 1) to have full rank 
conditioning on 74(A) having full rank. Again conditioning on 
A' having full rank, 



uo(fe-l) ui(fe-l) 

vo(^^-l) v,(^^-l) 
vo(/c' + A) vi(/c' + A) 



has probability at least 1 - 2{D + l)/(2^ - 1) to have full 
rank. Therefore A(A + 1) has full rank (i.e., equations in 
(28)U{zo(A; + A), ziik' + A)} are linearly independent) with 
probability at least (1 — 2(I^ + l)/(2^ — 1))^ conditioning on 
A (A) having full rank. 

Except for {2^0(^0), ^i(^i)}to=/c-i,ti=/c'-i, the equations 
of which are linearly independent with probability at least 
1 — 1/(2* — 1), all pairs of equations has probability at least 
(1 - 2(i:> + l)/(2* - 1))^ to be linearly independent with the 
previous pairs of equations. Therefore equations inside (21) 
are linearly independent with probability at least 

(,_^)(l_2(D + l) 



With the above lemmas, we can show that the broadcast 
throughput upper bound is achievable in grid networks with 
high probability. This result is presented in the following 
theorem. 

Theorem 2. The broadcast throughput in grid networks 
reaches 2/3 with high probability when s is of order larger 
than logD'^. 

Proof of Theorem 2: By Lemma 5, equations inside (21) 
are linearly independent with probability at least 



(1 



1 2(D+l),2^_2 



2s-i'^ 2«-l ' 
For large s, (1 — ^^^-i ^ ^2,0-2 ^.^^ approximated by 



(39) 



(40) 



Thus, if s is of order larger than logZ^^ (e.g., s = logD'^^), 
where e > 1, the limit of the above probability as ^ oo is 

1 



lim (1 



2(^+1) 2Z,_2 



1. 



(41) 



1^' 2^-1 

Therefore if D is large, at the end of round {D + 
max{/c,MA^ — k} — 2), node k can derive all native pack- 
ets from X with a high probability. At the end of round 
D + MN — 2, all nodes can derive all native packets from 
X with a high probability. The throughput is 

2D 2 
lim = — . 

D^oc 3{D + MN - 2) 3 



P 



(42) 



As a consequence, the broadcast throughput upper bound 
is achievable with high probability when the field size is of 
order larger than the logarithm of the number of packets. 

V. Conclusions 

In this work, we investigated the relationship between the 
broadcast throughput upper bound and min-cut in half-duplex 
broadcast networks. We showed that the theoretical throughput 
upper bound is n/(n + 1) for single-source broadcast, where 
n is the minimum vertex-cut size of the network. This upper 
bound is not always achievable in general, but is achievable 
in grid networks with the help of PNC. A ternary transmis- 
sion schedule with a three color coloring scheme has been 
proposed. By constructing a Hamiltonian cycle in the grid 
network, coloring along the Hamiltonian cycle, and schedul- 
ing nodes' transmissions based on the colors, the broadcast 
throughput is achievable when the field size is of order larger 
than the logarithm of the number of packets. In the future, we 
will extend our result and study the achievability of the upper 
bound in other networks. 

Appendix 

In this appendix, we argue that the broadcast throughput 
in the network in Fig. 2 cannot reach the upper bound given 
by Theorem 1. Recall that Wd is the number of time slots 
(or time window) before all nodes acquire all the D native 
packets. Let Ti be the set of time slots during which node i 
transmits within the Wd time slots; let 1Zi be the set of time 
slots during which node i receives within the Wd time slots. 

In order to achieve the the throughput upper bound (i.e., 
p = 2/3), we need 

m > ^iyz,,VzG{Xo,Xi, 1,2,3}. 
Since a node cannot transmit and receive at the same time. 



WD = \%unA = \n 



MTi\ =Wd- |7^.| < -Wd^i e {Xo,Xi,l,2,3}. (43) 
Consider node 1, the throughput for it is 

Pi < ITxo U Tx, I = ITxo I + ITx, I - ITxo n Tx, \ 

<lwD-\rxonTxA 



In order that pi > ^Wd, we need 

ITxonTxJ =0. (44) 

By focusing on node 2 and node 3, we can apply the same 
argument to show that 

\rxonT3\=0 (45) 

\Tx,nT2\=0 (46) 
The throughput of node 2 is 

P2<\rxo\^\r3\-\rxonr2\ m 

< lwD-\TxonT2\ (48) 
=>\TxonT2\=0 (49) 

if throughput of p2 is to be achieved. (47) is from the half- 
duplex constraint that when node 2 transmits during the T2 
slots, if Xq transmits at the same time, no information can be 
received. (48) is derived from (43). Similarly, we can argue 
that 

ITx, nTsl^O (50) 

|r2nr3| = o (51) 

Now, as Wd is the total number of time slots under consid- 
eration, we have 

^\Txo U Tsl + \Tx, U T2\ - \{Txo U T3) n {Tx, U T2)\ < Wd 
In order that p2 = ^Wd and ps = we need 

\rxours\>lwD 

\Tx,uT2\ > \wd 

The above are due to the fact that a node can receive 
information from the transmissions of its neighbors. 

l(rxour3)n(rx, ur2)| > \wd 

M{Txo n TxJ u (Txo n T2) u (Ta n TxJ u (Ta u T2)\ > ^Wd 
^\Txo n Tx, I + \Txo n Tsl + ITs n Tx, \ + ITs U Tsl > ^II^d 

^irxonrsi + iranrxj + iTsursi > ^ly^ (52) 

(52) is derived from (44). However, (52) contradicts (49), (50) 
and (51). Therefore, the throughput upper bound cannot be 
achieved. 
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